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PACS 05.70.Fh - Phase transitions: general studies. 
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Abstract - Inverse phase transitions are striking phenomena in which an apparently more or- 
dered state disorders under cooling. This behavior can naturally emerge in tricritical systems on 
heterogeneous networks and it is strongly enhanced by the presence of disassortative degree corre- 
lations. We show it both analytically and numerically, providing also a microscopic interpretation 
of inverse transitions in terms of freezing of sparse subgraphs and coupling renormalization. 
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Introduction. — Network-based representations are 
appropriate tools to describe many real- world systems and 
to classify their structural and functional properties using 
few topological measures, such as for instance the degree 
distribution, degree correlations and clustering PH. Net- 
works characterized by degree distributions with power- 
aw tails have attracted a lot of attention, mostly because 
of their unusual structural and dynamical properties [2jj4] . 
However these studies are mostly limited to uncorrelated 
random graphs, i.e. graphs characterized solely by their 
degree distribution. Many real-world networks show in- 
stead degree correlations. For instance, social networks 
are usually characterized by positive degree correlations 
(assortative mixing), because high-degree nodes tend to 
be directly connected and to form cliques [5]. In tech- 
nological and biological networks, instead, the hubs are 
preferentially connected to low degree nodes (disassorta- 
tive mixing) [6j. Degree correlations can strongly affect 
the overall behavior of dynamical processes, as it was re- 
cently shown for percolation [7] and diffusion processes 

In this Letter, we show that some amount of degree het- 
erogeneity and negative degree correlations can be respon- 
sible for the occurrence of inverse phase transitions. The 
latter are counterintuitive phenomena in which a system, 
starting from a high temperature disordered phase first un- 
dergoes a phase transition to a more ordered phase than 
comes back to a disordered one as a result of monotonic de- 
crease of temperature. Examples of this curious behavior 



are observed in a variety of systems, such as liquid binary 
mixtures, 3 He- 4 He isotopes, ultra-thin films, liquid crys- 
tals, disordered high-T superconductors and polymeric so- 
lutions [9j. An inverse transition implies the inversion of 
the standard ratio between the entropic content of the two 
phases 10 . This behavior has a simple microscopic expla- 



nation for some systems, such as water solutions of methyl- 
cellulose polymers, in which the more interacting unfolded 
state is entropically favored because it admits many more 
microscopic configurations than the non-interacting folded 
one. As suggested by the Flory- Huggins theory of poly- 



mer melts 11 , a general way of triggering an inverse phase 



transition is to introduce a temperature-dependent inter- 
action. This latter naturally emerges in spin systems with 



higher degeneracy of interacting states 12 . Moreover, an 



inverse freezing transition between glassy and paramag- 



netic phases exists in tricritical spin glasses 13 



and in spin-glasses on small world graphs 16 



14 



15 



Here we study the tricritical Ising model defined on 
random graphs, providing evidence of a novel topologi- 
cal mechanism that is responsible of an inverse melting 
transition. Using the cavity method, we show that degree 
fluctuations can generate a reentrance in the low temper- 
ature region of the phase diagram. The physical reason at 
the origin of this behavior is identified in the partial freez- 
ing of the microscopic degrees of freedom that generates 
a temperature-dependent renormalization of the effective 
interaction couplings. The mechanism, that requires some 
amount of degree heterogeneity to occur, is amplified by 
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Fig. 1: Phase diagram (r, 8) of the Blume-Capel model in the 
Curie- Weiss approximation for random graphs. The second 
order continuous branch has the universal scaling form of the 
Eqj3] The first order discontinuous branches are different for 
homogeneous (random regular graph of degree K — 5) and 
heterogeneous graphs (random graph with power-law degree 
distribution P(k) oc k~~* , with 7 = 3.5) respectively. 



disassortative mixing. This is of great relevance because 
many real-world technological and biological networks are 
disassortative and degree anti-correlations have been re- 
cently shown to be a natural feature of heterogeneous net- 
works (17) . 

The Model. — The tricritical Ising model, known as 
Blume-Capel (BC) model 



18 , describes spin-1 variables 



Si = {0, ±1} defined on the sites of a lattice or a graph, 
and interacting according to the Hamiltonian 



H B c = - j Y1 SiS i + A H ' 



(1) 



where the first sum runs over nearest neighbor pairs 
and the chemical potential A controls the density. It is 
well known that the BC model exhibits a rich T — A 
phase diagram with lines of first-order and second-order 
phase transitions dividing the ferromagnetic phase at low 
temperature T and chemical potential A from the param- 
agnetic phase at large T and/or A |l8j. The composition 
of the paramagnetic phase changes gradually and contin- 
uously from a purely disordered one (spins ±1) at large 
T and low A, to a vacancy-dominated state for large A 
and low T. In the following, we consider spins defined 
on the nodes of an uncorrelated random graph with de- 
gree sequence fci, &2, . . . , k^ randomly drawn from a given 
degree probability distribution P(k). 

Role of Degree Heterogeneity. — A qualitative in- 
sight on the role played by degree heterogeneity in deter- 
mining the critical behavior of the BC model can be ob- 
tained within a Curie- Weiss (CW) approximation. This is 
comparable to an "annealed" or a "random neighbor" ap- 



proximation where the adjacency matrix a^jrj is replaced 
in the Hamiltonian by the probability kikj/(k)N that the 
link between nodes i and j is present (2j. This approxi- 
mation allows the factorization of the free energy and the 
analytic computation of the partition function by Gaus- 
sian integration. It follows that the magnetization ^ of a 
neighbor of a randomly chosen site satisfies the mean-field 
self-consistent equation fi = F(jj) where 



k=k„ 



kP(k) sinh(^fc^) 
(k) e' 3A /2 + cosh(/3fc/x)' 



(2) 



This equation has a paramagnetic solution fi = whose 
stability is determined imposing F'(0) < 1. Following 
Landau theory, the equation F'(0) = 1 defines a curve of 
second-order critical points between the paramagnetic and 
the ferromagnetic phases (the A-line) till, upon decreasing 
the temperature, the third-order term in the free energjj^] 
expansion vanishes F"'(0) — (the tricritical point), then 
the transition becomes first order and the aforementioned 
condition gives us only the location of the spinodal point 
beyond which the paramagnetic solution is not stable. In 
terms of the rescaled variables 8 = A(k)/(k 2 ) and r = 
T(k)/(k 2 ) the A-line has the form 



5 = rlog[2(l/r-l)], 



(3) 



whereas the tricritical point is r c = | , 8 C = , that coin- 
cides with the maximum of the curve Eq(3j Interestingly, 
the obtained mean-field A-line depends on the graph prop- 
erties only through the ratio (k 2 )/(k), which determines 
the scaling factor. 

Lowering the temperature beyond the tricritical point 
(t c ,8 c ), we find a different dependence of the (first-order) 
transition line on the topological properties. It is useful 
to look directly to the T = limit, where Eq. [2] becomes 



H = {kQ(k/i-A))/{k) 



(4) 



and 0(x) is the Heaviside step function. A direct cal- 
culation, done approximating sums on k with integrals, 
reveals the following behavior. For graphs with a degree 
distribution which falls off faster than k~ 3 for large k, we 
find that the ferromagnetic state with \i = 1 is possible 
at T = when A < k m i n . In order to access the sta- 
bility property of this solution one needs to compare its 
energy E\p] — — (fc)/i 2 /2 + Afi with the paramagnetic en- 
ergy (E = 0) of the state with S{ = Vi. This shows 
that the fj, = 1 state is stable as long as A < (k)/2. So 
for graphs with k m i„ > (k)/2, as for example graphs with 
P(k) ~ fc~ 7 and 7 > 3, the coexistence region at T = 
extends from < A < k min with a first order phase tran- 
sition at (k)/2. For scale free graphs with 2 < 7 < 3 
the same arguments show that, at T = 0, a ferromag- 
netic state with /i ~ A 1- 2 exists for A < k max and it is 



Oj is 1 if i and j are neighbors, otherwise. 
2 This can be obtained by integrating eq[2] 
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Fig. 2: Reentrant behavior in a random graph of size N — 10 4 
and power-law degree distribution P(k) oc fc -7 . Left) Curves 
m(T) for different values of 7 and A. Right) Phase diagram 
for 7 = 2.5, 3, 3.5 obtained from the numerical solutions of the 
BP equations on single instances (points) and at the ensemble 
level (lines). 



thermodynamically stable for A < ((fc)/2) 7-2 . The same 
holds for 7 < 2 with the important difference that (k) now 
diverges with k max and that the magnetization fi is only 
weakly dependent on A. 

Some of the results are summarized in Fig (T] which dis- 
plays the phase diagram obtained in the CW approxi- 
mation for a homogeneous graph (regular random graph 
with degree K — 5) and for a heterogeneous graph with 
P(k) oc k~ 3 5 . The first-order branch was determined by 
equating the free energy of the paramagnetic and the fer- 
romagnetic solutions. While the lines of continuous phase 
transitions collapse using the rescaled variables (r, 5), the 
behaviors of the discontinuous transitions remain differ- 
ent as they depends on the average degree (fc), and not 
on the ratio (fc 2 ) / (k). The reentrant behavior of the first- 
order transition line is evident in the case of heterogeneous 
graphs. A similar result holds for queuing models on ran- 
dom graphs, where the transition point that divides con- 
gested and free phases scales differently with the moments 
of the degree distribution for continuous and discontinu- 
ous transitions, respectively |19| . In summary, the CW 
approximation would suggest that in the BC model an 
inverse phase transition can be triggered by a different 
scaling of first and second order critical lines with the mo- 
ments of the degree distribution. We see in the following 
that, although the role of degree heterogeneity is crucial 
to have a reentrance, this is not strictly associated to the 
different scaling of the first-order and second-order phase 
transitions. 



Results in the BP Approximation. The exis- 
tence of a topology induced inverse behavior is confirmed 
by the more accurate Bethe-Peierls (BP) approximation 
[2]. On a tree-like structure emerging from a randomly 
chosen node i, the partition function of the system can be 
written as 



n m *j( s *) 



(5) 



where &i is the set of neighbors of i and rriij(si) satisfies 
the recursive equation 



mij(si) 



(6) 



k£dj\i 



we 



With the parametrization mij(si) = Aye^ Ui 
can easily obtain a set of recursive equations for the fields 

m 

s 2 ) and local 



Uij and the anisotropics Vij , that can be solved numerically 
and used to compute the local densities p 
magnetization mi — (sj) for each node i: 

2cosh(/3£ 



Pi 



rn, 



3 U V' 



e /3(A+2: 3 v tj ) _ 
Pi tanh(^ u t 



2cosh(/?£. 
))■ 



(7) 



In Fig{2] (left) we plot the total magnetization as a 
function of the temperature T in random graphs with 
P(k) oc fc~ 7 (N — 10 4 ). For all values of 7, the magnetiza- 
tion at T = drops to zero for sufficiently large A, even if 
a metastable ferromagnetic phase survives at T = 0. The 
agreement between BP results (lines) and Monte Carlo 
simulations (points) is very good. In the right panel of 
Fig[2j we report the phase diagram obtained using the BP 
approximation for the same values of 7. The points are 
obtained averaging the results obtained running BP equa- 
tions on given instances, while the full line is the analytic 
prediction of the continuous transitions and the tricriti- 
cal points obtained again within the BP approximation, 
but solving the equations at the level of ensembles of ran- 
dom graphs. Indeed, for random graphs characterized only 
by the degree distribution P(k) we can divide the nodes 
in classes according to their degree k and consider the 
cavity fields symmetric within each class. Once 

introduced the average cavity fields u — ^ fc P{k)uk and 
v — J^k P(k)vk, one can easily check that they have to sat- 
isfy a closed set of self-consistent equations, u — f(u,v), 
v — g(u, v). The paramagnetic solution of these equations 
reads 



u = 



e /3(A+(fe-l)5*) + 2 
,f}(A+(k-l)v*) +2 cosh ]) 



(8) 



Upon expansion around this solution, the tricritical point 
is given by the conditions d n f(0,v*) = 1, 8^/(0, v*) = 0. 
The first equation describes the A line of the second order 
critical points and gives us the value of A c (/3) as a function 
of the temperature while the second equation fixes the 
ending point of this line, /3 C and its associated A C (/3 C ). 
The previous conditions result into a set of two coupled 
implicit equations for (3 C and A c 

fc(fc-l)P(fc) 2sinh/3 c 



E 



(k) 



E 



k{k - lfP(k) e 



2f,. 



e l " + 2 cosh / 
2e ia cosh fl c — 



(k) 



{e l - +2cosh/3 c 



1 



(9) 
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Fig. 3: The free energy computed on a single graph is shown 
for A = 4 in the left panel and for A = 3 in right panel. We use 
the algorithm developed in 25 , by using a Legendre transfor- 



mation. The scale free graph used has N = 5000, k min = 4 and 
7 = 3.5. In both cases the free energy displays a minimum in 
zero for T > T c the critical temperature and, more importantly, 
for T < T m (the melting temperature). In between these two 
temperatures, T m < T < T c , it shows a ferromagnetic behavior 
with two minima in ±m. 



where e 



/3 c (A c + (k-l)v') 



These equations have been 



solved numerically for ensembles of power-law random 
graphs with 7 = 2.5, 3, 3.5 (see right panel in Fig. |2|. 

The overall picture obtained from the CW approxima- 
tion is confirmed by BP results, although the position of 
the tricritical point seems to move to the left of the maxi- 
mum of the curve in the BP phase diagram, showing that 
inverse transitions can be continuous as well. 

We checked that the fixed point solutions of the cavity 
equations have a true thermodynamic meaning by com- 
puting and plotting the Landau free energy as a function 
of the magnetization f(m). We use the algorithm devel- 
oped in [25] to reconstruct the free energy landscape by 
means of a Legendre transformation. The two panels re- 
ported in Fig|3] show that the free energy of a scale free 
graph for A = 4 (left) and A = 3 (right). We can see 
that in both cases the free energy signals the presence of a 
stable paramagnetic solution at high temperature, i.e. it 
has a minimum at m = 0. Then a ferromagnetic solution 
m > appears at smaller temperature and finally, cooling 
down the system below the melting temperature T m , it is 
possible to see the minimum of free energy again at the 
origin m « 0. 

A Microscopic Mechanism for Reentrance. By 

performing a low temperature analysis, we found that 
nodes with different degrees can have very different val- 
ues of magnetization. In FigQtop-right) we report the 
magnetization curves m(T) at A = 3 for a bimodal ran- 
dom graph with P(k) = 0.2^^0 + 0.85k,2 obtained by 
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Fig. 4: Left) Effective interaction J e ff(T) between the neigh- 
bors of a node of degree 2 as a function of the temperature T for 
several A. Right) Magnetization as a function of the temper- 
ature for A = 3 in uncorrelated (top) and disassortative (bot- 
tom) bi- modal random graphs with P(k) = 0.2(5^ 2 + 0.85fc, 10. 
The behavior of nodes belonging to different classes of degrees 
is highlighted. The curves display results from both Monte 
Carlo simulations (points) and BP calculations (solid lines). 



Monte Carlo simulations and BP calculations. The curve 
referring to nodes with lower degree is showing a reen- 
trant behavior, whereas the one of high degree nodes is 
not. This behavior suggests that the topology- induced 
reentrance could be activated at low T by a complex re- 
cursive microscopic mechanism starting from the spins of 
lowest degree that freeze to zero for a sufficiently large 
chemical potential A. For simplicity we consider a node 
i of degree 2 and two neighboring nodes ii, 12 of degree 
larger than 2. At zero-temperature a chemical potential 
A > 2 is sufficient to freeze s, = independently of the 
neighbors of i. Once Sj is frozen to zero, the two neighbors 
i±, 12 are not directly correlated anymore, the effective in- 
teraction is zero and the local fields acting on them do 
decrease. They now have to be compared with A to es- 
tablish whether nodes i\, 12 are set to zero or not. Pushing 
forward this decimation procedure (which coincides with 
A-core percolation [20]) makes it possible to recursively 
determine the whole ground state of the system. We can 
also check how the interaction between Si 1 and Si 2 changes 
at finite T. If we sum over the values of the internal spin 
Si, the partial partition function of the three spins be- 
comes that of a system of two spins with an effective cou- 



pling J eff = ± log f 1+2e 1+ P 2e -° S A (2/3) j ■ In the left part 
of Fig. [4] we show that as soon as A > 1, the coupling 
between i\ and 12 becomes non-monotonic in T and for 
A > 2 the coupling vanishes at T = due to the freezing 
of Si = 0. Combining the renormalization argument with 
the previous decimation algorithm, it is possible to check 
that two far-away spins remain completely uncorrelated 
at low T and their effective coupling increases at larger 
temperatures. Degree heterogeneity is the first essential 
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ingredient needed for an inverse transition because when 
all nodes have approximately the same degree, the renor- 
malization of couplings is not effective. Another impor- 
tant ingredient that enhances the reentrant behavior are 
degree anticorrelations. In order to verify this idea, we 



used a Monte Carlo algorithm proposed in 22 in order to 



generate instances of correlated networks where high de- 
gree nodes are preferentially connected to low degree ones. 
An uncorrelated random graph is first generated with the 
configuration model. Then Monte Carlo moves are used to 
exchange the ends of randomly chosen link-pairs, in order 
to introduce positive or negative correlations among the 
degrees without altering the degree distribution. In the 
bottom-right panel of Figj4]we show the curves m(T) for 
the BC model in a bimodal random graph with disassorta- 
tive mixing. In this case the freezing of low degree nodes 
at low T is sufficient to disconnect the whole graph trig- 
gering an inverse phase transition. Similarly, a low level of 
disassortative mixing is sufficient to enhance the reentrant 
behavior in scale- free graphs (not shown). Therefore, dis- 
assortative mixing favors inverse phase transitions. Unlike 
mean-field results, this argument is not affected by spatial 
fluctuations, hence we guess that a topologically induced 
inverse transition should be observable in dimensional lat- 
tices as well if there is a sufficient degree of local disorder 
and degree anti-correlation. The renormalization argu- 
ment as well as the mean-field calculation can be extended 
to other tricritical spin models, such as the random-field 
Ising model (RFIM) with a bimodal distribution of fields. 
However, the phase diagram of this model has a complex 
structure of singularities even on homogeneous random 



the generality of the observed phenomenon. 

Conclusions. — We showed that tricritical spin sys- 
tems can admit inverse phase transitions on heterogenous 
graphs. Mean-field arguments qualitatively suggest a dif- 
ferent scaling of the first and second order critical point 
with the moments of the degree distribution for random 
graphs. A microscopic mechanism of coupling renormal- 
ization shows that in the presence of degree heterogeneity, 
the effective coupling among two nodes becomes temper- 
ature dependent. In this way, we have shown that the 
conditions usually described in the Flory-Huggins theory 
of inverse melting can have also a topological nature. The 
reentrant behavior is strongly enhanced (damped) by the 
presence of negative (positive) degree correlations. 



DDM warmly thanks Luca Leuzzi for fruitful discussions 
and useful suggestions. 
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